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ABSTRACT 


Several mechanistic models are developed to predict recombination coefficients for use in 
heat shield design for reusable surface insulation (RSI) on aerobraking vehicles such as 
space shuttles. The models are applied over a temperature range of 300 to 1800 K and 
a stagnation pressure range of 0 to 3,000 Pa. A four parameter model in temperature 
was found to work best, however several models (including those with atom concentrations 
at the surface) were also investigated. Mechanistic models developed with atom 
concentration terms may be applicable when sufficient data becomes available. This work 
shows the requirement for recombination experiments in the 300 to 1000 K and 1500 to 
1850 K temperature range, with deliberate concentration variations. 
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INTRODUCTION 


The objective was to determine a continuous model for the recombination coefficients of 
atoms on space shuttle tiles and similar heat shield materials. This model could then be 
applied in heat shield design calculations used in reentry simulations, to help determine 
the amount and type of heat shield material required. Currently, the experimental data 
shows a maximum recombination rate at around 1600 K, accounted for by discontinuous 
recombination models which change at the maximum observed in the data. This work 
closely examines two mechanisms to determine if they can predict the observed 
recombination data over a wider range of temperatures. These mechanisms also 
demonstrate an influence of atom concentration on the recombination coefficients. 

The results can be applied in heat shield calculations for the aeroassist flight experiment 
(AFE) and aero-orbiter transfer vehicle (AOTV). Higher reentry temperatures are 
predicted at the heat shield surfaces. However, current recombination models cannot 
predict observed recombination rates at higher temperatures. For instance, a simple 
Arrhenius model predicts heating rates two times higher than measured because the 
model predicts higher recombination coefficients than those observed. The use of a 
continuous model'which reflects actual recombination data across a wider temperature 
range, will allow more flexibility in determining the best reentry pattern and heat shield 
design for space craft reentiy systems. 

Figure 1 from Ting et al. (1989) demonstrates the aeroassist flight experiment (AFE) 
which is proposed for 1993. This figure portrays how the aeroassist vehicle will travel 
through the atmosphere, simulating an aerobraking maneuver while returning to the space 
shuttle. It is during this aerobraking maneuver where recombination and heat shield 
design become important. 
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BACKGROUND 


General Background 

Space shuttle orbiter reentry experiments have shown that catalytic surface recombination 
influence heat shield temperatures (Scott, 1985). Specifically, increases up to 111 K have 
been reported for temperatures measured on a catalytically coated space shuttle tile in 
contrast to non-catalytically coated space shuttle tiles. The higher temperatures observed 
on catalytically coated tiles are explained by increased surface recombination of atoms to 
molecules. Further, Rosner & Fang (1974) suggested that recombination is not entirely 
to ground state molecules. For instance, molecules created by surface recombination at 
one point come off the surface excited, later transferring energy further downside the 
space craft surface. To date, mechanistic description of this phenomena has been lacking. 
Based on data presently available, a mechanistic model should be able to describe the 
maximum in the recombination coefficient as temperature increases. This maximum occurs 
around 1600 K based on available data. 

Figure 2 demonstrates the reentry chemistry that occurs. As a vehicle reenters the 
atmosphere, a bow shock is established off the front nose of the vehicle. As flow goes 
across the bow shock, molecules are disassociated into atoms. These atoms then 
recombine on the surface, dumping their energy of recombination onto the surface, and 
leave the surface as molecules. The ideal surface would have four characteristics: 1) high 
emissivity, 2) low recombination rates, 3) the ability to withstand high temperatures, and 
4) insulating capabilities to the vehicle itself. This report focuses on the description of 
recombination rates of atoms. To demonstrate the amount of energy released during 
recombination, it is important to note that the reaction of two oxygen atoms on the 
surface will release about 500 kJ/mole. This is equivalent to the energy released by the 
combustion of coal on a per unit mass basis. Similarly, recombination of nitrogen atoms 
releases about two times as much energy on the surface. Therefore, each pound of 
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FIGURE 2 


Reentry Flow Field 
(from Kolodziej and Stewart, 1987) 
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nitrogen atoms which recombine on the surface would be equivalent to burning two 
pounds of coal. Space shuttle tiles are composed of about 1" thick low density silica called 
"high temperature re-usable surface insulation" (HRSI) coated with about 3 millimeters 
thick reaction cured glass (RCG). RCG coating is an excellent non-catalytic surface and 
is ideal for the purpose of low recombination. However, some recombination continues 
to occur and must be accounted for in heat shield calculations. 

Currently, recombination is accounted for by recombination coefficients which are the 
combination of two components. The first component is the probability that atoms will 
recombine when they collide on the surface. The second component, then, is the fraction 
of dissociation energy released to the surface. The recombination coefficients discussed 
in this report are thus called energy transfer catalytic recombination coefficients for 
nitrogen and oxygen. The ability to separate the probability of recombination from the 
fraction of dissociation energy released has yet to have been accomplished experimentally. 
However, mechanisms can still be used to describe the data. In general, recombination 
data is fitted to an Arrhenius expression: 

* = k (1) 

where k consists of 2 parameters, k 0 and Ea. The Arrhenius expression becomes: 
k = k 0 exp (-Ea/RT) (-la) 

where k 0 is called the pre-exponential coefficient and Ea is called the activation energy. 

For example, Kolodziej & Stewart (1987) reported the recombination coefficient for 
nitrogen and oxygen as shown below. The Arrhenius expression changed at about 1600 
K for nitrogen and oxygen, and specifically the activation energy changed from positive to 
negative. 


N = 6.1 x 10' 2 exp (-2480AV) 

1410 <T W < 1640 

0 

* „ = 6.1 x 10‘ 4 exp (5090/TV) 

1640 <T W < 1905 

(3) 

Y 0 = 40 exp (-11440/T^) 

1435 <T^< 1580 

(4) 

- 39 x 10' 9 exp (21410/T^) 

1580 <Tjp< 1845 

© 
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Recombination data taken over a wider range of temperature suggest that a simple 
Arrhenius expression will not work. For example, look at the recombination coefficient 
of nitrogen on silica in Figure 3 (Newman, 1987). Ignoring the wide range of data at 300 
K, the figure shows that from a temperature range of about 300 to 1000 K a rather low 
activation energy exists (proportional to the negative slope of the line). As temperature 
increases from 1000 to about 1400 K the activation energy increases. Above 1400 K, the 
activation energy changes from positive to negative and gamma decreases. 

The objective of this work is to describe the recombination coefficient over a wider 
temperature range with one equation. 

Recombination Data Background 

Recombination of atoms on RCG coated HRSI surfaces was the data used in this 
mechanistic modelling study. The data came from six sources: Kolodziej and Stewart 
(1987); Scott (1981); Stewart (oxygen data only) (1982); Willey (summer 1988) and 
Marinelli (1986) & (1988). Kolodziej, Scott, Stewart & Willey acquired their data from 
arc plasma test facilities. In Kolodziej’s work, four RCG coated HRSI models, with a 
thickness of 3 mm., were used. The HRSI material was FRCI 20. Temperature 
measurements were made with Pt/PT-13% Rh thermocouples and the emissivity was 
reported at 0.89. Heating rates were determined by a calorimeter of known catalytic 
activity. To predict the recombination coefficient, the heating rate of a known catalytic 
reference was measured compared to the heating rate as determined by the Stephan- 
Boltzman equation, Q = e<TTV The calculation method used by Kolodziej and Stewart 
(1987), a simulation modeling and analysis of reusable thermal protection package 
(SMART), consisted of 3 programs: NOZZL, which calculates properties of a hot gas 
expanding through an arc jet nozzle; AMIR which models the inviscid bowshock in front 
of the blunt body; and BLIMP which models the non-equilibrium boundary layer with 
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FIGURE 3 


Recombination Coefficients for Atomic Nitrogen 
on Silica Surfaces (from Newman, 1987) 

Temperature, K 
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surface reaction kinetics. These programs generate curves of Q observed/Q fully catalytic 
(Q/Q FC ) as a function of gamma. Gamma is determined by the intersection of the 
experimental data with the predicted curves. More details are described in Kolodziej et 
al. 1987. 

Scott and Willey also used arc plasma facilities to measure recombination coefficients. 
Scott (1981) used Pt/Pt-10% Rh thermocouples with the gas compositions of nitrogen and 
air; the enthalpies ranged from 20-37 MJ/kg and 5-31 MJ/kg for nitrogen and air 
respectively. The model was a flat faced RCG coated HRSI model. Willey (1988) used 
conditions similar to Scott with the exception that oxygen compositions were varied from 
0, 11, 14 & 21 mole percents. A hemisphere model with a thermocouple located at the 
stagnation point of the hemisphere was used. The calculational approach used by both 
Scott and Willey is also described in Scott (1981). Flow field properties are determined 
by simple flow conservation equations for an inviscid stream tube and the flow is assumed 
to be a chemically frozen ideal gas. Q/Q FC versus gamma curves are generated and an 
experimental gamma is determined from these curves. More details about the 
calculational routine used by Scott and Willey is available in Scott (1981). 

The work of Marinelli was undertaken in a discharge flow reactor. Atom concentrations 
were determined by using laser-induced fluorescence and resonantly enhanced multiphoton 
ionization spectroscopy. Gamma was determined by solving the first order kinetic 
equations for oxygen concentration and was proportional to light emitted at 225 nm. More 
details about Marinelli experiments can be found in Marinelli, 1988. 
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DATA USED IN THE ANALYSIS 


The. mechanistic models developed in this work, were fitted to recombination coefficient 
data shown in Tables 1 and 2, with the first column representing the source of the data. 
Kolodziej and Stewart’s data are from their 1987 AIAA paper. The data from Scott were 
presented in his 1981 paper. Two points for nitrogen from Scott were deleted from this 
analysis because they were off the trend shown by the bulk of his data. The data of 
Willey were acquired during the summer faculty fellowship program at NASA-JSC during 
the summer of 1988 and was reported in a summer faculty fellowship report. The datum 
for Marinelli (1986) for nitrogen atoms was found in a secondary source in a thesis by 
Newman, 1987. The second column in Tables 1 & 2 presents the temperature as read 
from the figures in the paper or reported. The third column reports the estimated partial 
pressure of the atoms at the surface of the RCG coated test article. Partial pressure was 
equal to the total stagnation pressure times the atomic mole fraction as reported or 
estimated in the various sources. The fourth column demonstrates the recombination 
coefficient (gamma) measured or reported. Again, gamma represents the energy transfer 
recombination coefficient which is a combination of the probability that atoms will combine 
times the fraction of energy released. 

Atom concentration is related to the partial pressure by the ideal gas law: [O] = P 0 /RT. 
In the reaction models which follow, the concentration of oxygen was used in place of the 
partial pressure of atoms. 

Figures 4 and 5 show the log of gamma plotted against 1/T for all data with the exception 
of Marinelli’s data. The data show the general trend of reaching a maximum at about 
1600 K, (1000/T=0.667). Gamma decreases as temperature increases above a temperature 
of 1600 K (looking from right to left in the figures). 
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TABLE 1 


Raw Data Used in the Estimation of Parameters for Various 
Mechanistic Models for the Recombination of Nitrogen Atoms on 
RCG Coated HRSI Surfaces 


Source 


Kolodziej & Stewart 
(1987) 


Scott (1981) 


Willey 


(1988) 


Marinelli (1986) 


Temp., 

K 

Partial 

Pressure 

Nitrogen 

Pa 


1889 

2073 

0.00894 

1742 

2073 

0.00943 

1686 

608 

0.01380 

1603 

608 

0.01465 

1467 

608 

0.01194 

1422 

608 

0.01068 

1668 

1010 

0.01945 

1661 

840 

0.01478 

1582 

721 

0.01740 

1546 

751 

0.01499 

1438 

461 

0.02155 

1387 

521 

0.01534 

1092 

170 

0.00943 

1022 

128 

0.00810 

1014 

120 

0.00772 

996 

92 

0.00687 

1536 

291 

0.01504 

1524 

295 

0.01401 

1506 

279 

0.01267 

1496 

257 

0.01723 

1493 

268 

0.01464 

1490 

266 

0.01527 

1460 

240 

0.01156 

1412 

205 

0.01383 

1411 

208 

0.01437 

1402 

180 

0.01368 

1350 

207 

0.01540 

300 

133 

0.00019 
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TABLE 2 


Raw Data Used in the Estimation of Parameters for Various 


Mechanistic Models for the Recombination of Oxygen Atoms on 


RCG Coated HRSI Surfaces 


Partial 

Temp Pressure 

Source K Oxygen 0 O 

Pa 


Kolodziej & Stewart 

1831 

824 

0.00370 

(1987) 

1806 

412 

0.01007 

II 

1742 

412 

0.01893 

ft 

1726 

824 

0.00525 

1? 

1644 

235 

0.01068 

fl 

1617 

412 

0.00968 

II 

1592 

235 

0.03060 

H 

1450 

235 

0.01790 

Scott (1981) 

1647 

370 

0.02361 

ii 

1493 

368 

0.01300 

ii 

1419 

351 

0.00775 

Stewart et al. (1982) 

1556 

451 

0.01733 

M 

1344 

406 

0.01120 

Willey (1988) 

1557 

241 

0.00887 

if 

1552 

134 

0.02793 

M 

1540 

252 

0.00885 

M 

1527 

253 

0.01081 

II 

1522 

127 

0.01637 

II 

1518 

241 

0.00817 

ii 

1501 

237 

0.01524 

II 

1497 

239 

0.00838 

ir 

1486 

230 

0.02148 

m 

1479 

128 

0.01972 

ii 

1458 

225 

0.02432 

n 

1456 

124 

0.01146 

li 

1449 

231 

0.01500 

Marinelli (1988) 

300 

133 

0.00020 
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RESULTS AND DISCUSSION 


Models Evaluated 

Arrhenius Expression 

The first model compared to the data was a simple Arrhenius expression given by: 

2fo = 2k (6) 

where k is a function of temperature given by: 

k = k 0 exp (-Ea/RT) (6a) 

and where k 0 = pre-exponential constant and 
Ea = activation energy. 

The value of 2 in Equation 6 reflects the stoichiometric requirement that two times the 
rate of molecular oxygen production equals the rate of atomic oxygen consumption. This 
model demonstrates that chemical reactions contain a barrier between reactants and 
products. This barrier is called an activation energy, or Ea. The higher the Ea the more 
difficult it is for reactants to go to products at low temperatures. The reaction rate is 
much faster at higher temperatures because higher temperatures provide the energy 
needed to hurdle the higher barrier. The Arrhenius model is a direct result from a 
standard definition of the recombination reaction being first order in oxygen atoms. The 
recombination rate is equal to k H , times the concentration [O]. 

rate = k*, [O] (7) 

The recombination coefficient is defined as the ratio of the number of atoms recombined, 
divided by the number of atoms that strike the surface. The number of atoms which 
strike the surface is given by: 

N = [O] c/4 = [O] [k B T/2m7f , ] ;/2 (8) 

where c is the average velocity per particle equal [8k B T/mT/] //2 
and [O] is the number density of atoms, atoms/m 3 . 
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Dividing Equation. 7 by Equation 8 gives the definition of gamma in terms of the first 
order wall reaction coefficient, k^. 

% = K(2m^ B T) m (8a) 

Thus, the constant k in Equation 6 is similar to the constant k*, in Equation 8a within a 
correction factor related to the square root of temperature. Often in chemical kinetics the 
square root of temperature relationships are assumed to be constant coihpared to the 
exponential factors in the reaction coefficients. Thus, the recombination coefficient for the 
Arrhenius expression assumes that the reaction rate at the surface is first order in oxygen 
concentration above the surface. 7 The Arrhenius expression predicts that when the log 
of gamma is plotted against 1 /T, a straight line will result, whose negative slope is 
proportional to the activation energy and whose intercept is related to the preexponential 
constant. Reviewing the data in Figures 4 and 5, one sees that, indeed, a portion of the 
data can be described by a straight line going from right to left of a negative slope. After 
the temperature reaches about 1600 K, however, the Arrhenius expression fails to describe 
the data properly. The Arrhenius expression will be referred to as Model 0. 

Recombination Models Derived from the Rideal-Elev Mechanism 
The Rideal-Eley mechanism is a two-step mechanism. First an atom adsorbs onto a 
surface site, then the adsorbed surface atom reacts at the surface with a gas phase atom. 
The mechanism is as follows: 

O + [S] [OS] (9) 

O + [OS] — > O, + [S] (10) 

In Equation 9 , oxygen atoms are adsorbed onto reduced sites. In Equation 10, oxygen 
atoms near the surface combine with oxygen atoms on the surface to form O,, plus 
regenerate a reduced site. Several reaction models can be found from this mechanism and 


1 For the recombination of nitrogen atoms, substitute nitrogen concentrations 
whenever oxygen concentration is discussed. 
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are derived in Appendix A. The simplest model from the Rideal-Eley mechanism, is to 
assume that Equation 9 is in pseudo-equilibrium and that Equation 10 is rate limiting. 
Thus, the simplest model that can be derived is as follows: 

= 2k 2 K/(l + K) (11) 

This equation will be referred to as Model l 2 and, when compared to Model 0, has an 
added term related to adsorption (K). This model predicts that at low temperatures, 
where adsorption dominates, gamma will be zero order in oxygen atom concentration. 
Equation 11 reduces to: 

0 = 2k, (12) 

The surface will be essentially bare above a certain temperature, demonstrating that as 
temperature increases, fewer oxygen atoms will stay on the surface. The reaction rate can 
then be approximated by: 

= 2 k, K. (13) 

The adsorption term, K, consists of the ratio of the forward divided by the reversed 
reaction rate constants in Equation 9. In general, the adsorption equilibrium constant is 
described by an Arrhenius behavior as given below: 

K = Ka exp (D/RT) (14) 

where Ka is the preexponential constant for the equilibrium constant and D is the 
activation energy related to the adsorption equilibrium constant. 

The equilibrium constant has a positive slope on an Arrhenius plot because adsorption is 
stronger at lower temperatures. As temperature increases, molecules tend to come off the 
surface at a faster rate than they can be absorbed. The overall results observed are 
dependent upon the value of D. The combined k 2 K term can have either a positive or 
negative overall activation energy. Therefore, the behavior of the recombination 
coefficient at higher temperatures, can increase, stay the same, or decrease in temperature 
depending on the value of D compared to Ea. 


2 Derivations of Models 1 through 4 are in Appendix A & B. 
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Model 1 assumes that the concentration of oxygen atoms near the surface is constant for 
all conditions. Another Rideal-Eley model, Model 2, includes a concentration of oxygen 
atoms at the surface and is as follows: 

% o = 2 k 2 K[0J/(l+K[0]) (15) 

This model shows that oxygen will influence the recombination rate. At low temperatures, 
Model 2 reduces to: 

£» = 2k 2 (16) 

and at high temperatures, Model 2 reduces to: 

= 2k,K[0] (17) 

Therefore, the recombination order will be zero order at low temperatures transitioning 
over to first order at high temperatures. 

Model 2 effectively displays atom concentrations at the surface proving that recombination 
coefficients are a function of atom concentrations near the surface. Unfortunately, there 
is not enough data currently available to determine the true functionality of recombination 
coefficients with surface concentrations. 

An extension of Model 2 includes a temperature correction term which appears in the 
derivation of the model: 

(2rntf[k D T) I/2 (IS) 

Model 3 will essentially have the same behaviors at the low and high temperature 
extremes as displayed by Model 2. Model 3 is as follows: 

= 2k 2 (2m7T/k B T) 7/2 K[0]/(l +K[0]) 2 (19) 
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A Recombination Model Derived from the Lanemuir-Hinshelwood Mechanism 
In the Langmuir-Hinshelwood mechanism oxygen atoms adsorb onto the surface to form 
a surface oxide species. Next, two neighboring surface oxide species combine to form 
molecular oxygen: 

O + [S] [OS] (20) 

[OS] + [OS] — > O, + 2 [S] (21) 

Equation 20 represents the adsorption of oxygen onto the surface and is similar to 
Equation 9 in the Rideal-Eley mechanism. Equation 21 represents the reaction of two 
surface oxide species to form molecular oxygen. The resultant recombination model, 
Model 4, is based on the assumption that Equation 20 is in pseudo-equilibrium and that 
Equation 21 is the rate limiting step. 

= 2k, (2mr/k fl T) ;/ -K-[0]/(l + K[0]) 2 (22) 

Model 4 demonstrates that at low temperatures, where adsorption dominates, the 
recombination coefficient is negative first order in oxygen concentration: 
fo = 2k 2 (2irfy/k s T)" 2 /[0] (23) 

At high temperatures, where adsorption is negligible, the recombination coefficient is first 
order in oxygen: 

y„ = 2 k,(2mir/k„T)" 2 K 2 [0] (24) 

Thus, as temperature increases, a transition occurs from negative first order to positive 
first order in oxygen atom concentration at the surface. The available data, however, are 
limited and more testing is suggested to verify the possibility of this mechanism. Once 
again, the constant, k„ and equilibrium constant, K, follow Arrhenius relationships as 
discussed in the Rideal-Eley mechanism. 
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Seward Model 

The Seward model was developed in his 1985 thesis and is based around the Rideal-Eley 
mechanism using first principle kinetics to describe the adsorption, desorption, and reaction 
rate steps separately. This model, Model 5, consists of eight parameters: 

^ 0 = 2PN 0 So exp(-Ea/RT)/(SoN 0 +g+ PN C exp (-Ea/RT)) (25) 

The term P represents the steric factor and is related to the ability of oxygen atoms to 
absorb. Seward allowed this factor as a function of temperature, up to a certain constant. 
The So term is related to the sticking coefficient and predicts the probability that an atom 
will adhere once it strikes the surface. Seward also demonstrated this constant as a 
function of temperature. The exponential term in the numerator is related to the forward 
reaction rate of step 2 in the Rideal-Eley mechanism. The term delta in Seward’s 
equation is related to the desorption step in Equation 9 of the Rideal-Eley mechanism and 
consists of an activation energy and several constants related to the number of surface 
sites. Further explanation of the Seward model is presented in Appendix A. 

Model 5 has three characteristic temperature ranges: at low temperatures adsorption of 
atoms onto the surface dominates, and reduces to: 

^ 0 = 2P exp (-Ea/RT) (26) 

As temperature increases, the forward reaction rate of step 2 increases. Thus, the 
exponential term dominates in the denominator and the numerator: 

X = 2So (27) 

Finally, at high temperatures, the desorption of oxygen atoms dominates and the following 
approximation results: 

^ = 2P N 0 So exp (-Ea/RT)/ g (28) 

Model 5 predicts zero reaction order at low and medium temperatures in oxygen atom 
concentration; at higher temperatures, the model becomes first order in oxygen atom 
concentrations. The steric factor and the sticking coefficient in Model 5 are allowed to 
be a function of temperatures, and therefore, at lower temperatures the model concaves 
upwards. An attractive feature of this model is that it predicts low activation energy for 
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recombination at low temperatures, changing to a higher activation energy at moderate 
temperatures. The difficulty in the application of Model 5 is the determination of the 
eight parameters (there are too many parameters for the data that is currently available). 

Willey Empirical Model 

A six parameter empirical model, Model 6, is compared to the data available: 

)S Q = 2(k 2 + kj)K/(l+K) (29) 

This model can be thought of as the combination of two mechanistic steps, the rate being 
determined by the fastest of the two steps, with each step dependent upon activation 
energy. In this work, the rate constant k 2 dominates at low temperatures. At moderate 
temperatures the recombination coefficient is described by k 2 and at high temperatures the 
recombination coefficient is described by 2k?K 

Model 6 can predict a concave upwards curve at low temperatures for the recombination 
coefficient as a function of 1/T. For the recombination of nitrogen atoms on silica, as 
shown in Figure 3, the activation energy for the recombination coefficient increases as 
temperature increases. Finally, at high temperatures, where desorption dominates, the 
overall activation energy goes from negative to positive and therefore the recombination 
coefficient declines as temperature further increases. Model 6’s parameters can be 
adjusted to reflect the behavior shown in Figure 3. Information concerning concentration 
dependence on recombination coefficients is lacking and, therefore, a concentration term 
is not included in Model 6. A summary of mechanistic models evaluated in this report 
is shown in Table 3. 
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TABLE 3 


Summary Mechanistic Models Evaluated 
for Recombination Coefficients for N and O Atoms. 

Model # Name Model 


0 

1 


2 


3 


4 


5 


Arrhenius ^f 0 = 2 k 2 

Rideal-Eley 2 k 2 K 

without ConcentrationV = 

Terms (1 + K ) 


Rideal-Eley 



2 k 2 K [O] 


( 1 + K [O] ) 


Rideal-Eley 2 k 2 (2mTj7k B T) ;/2 K [O] 

with a Temperature ^ = 

Correction Factor (+ K [O] ) 

Langmuir- 2 k 2 (2m7f7 k B T) m K- [O] 

Hinshelwood = 

( 1 + K [O] ) 2 


Seward 



2 P N a So exp(-Ea/RT) 

(SoN 0 + $ + P N 0 exp(-Ea/RT)) 


6 


Willey 



2 K (k, + k 5 ) 
( 1 + K) 


Notes: 

k 2 is a function of temperature given by: k 2 = k 20 exp(-Ea/RT) 


kj is a function of temperature given by: k 2 = k 30 exp(-Eaj/RT) 


K is a function of temperature given by: K = Ka exp(D/RT) 

[O] is the concentration of oxygen atoms at the surface in atoms/m 3 ; 
for nitrogen substitute [N] for [O] in any of the above models 


The value of R in these models is 0.0083144 kJ/g mole/K 
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Detennination of Constants 


For the models presented above, the parameters were determined by non-linear least 
squares on the log of gamma. The log of gamma was chosen because it would give equal 
weight to lower gammas at lower temperatures. Therefore, the parameters chosen were 
those which minimized the sum of squares of error (SSE) defined as: 2[(log^ w 0&s ) - log 
($ (i) pred)] 2 - All plots presented are log gamma versus 1/T. Model 0 was simple to 
evaluate by linear regression on log gamma versus 1/T because the slope is proportional 
to the activation energy and the intercept is equal to 1/2 of k 2 . The parameters for 
Models 1-4 were determined by using a non-linear least squares routine, BMDP PAR, for 
derivative free non-linear models. Models 1-4, it was discovered, have poorly conditioned 
surfaces and the parameter’s joint confidence regions have very narrow but broad cigar 
shape surfaces (Draper & Smith, 1986). To provide assistance in determination of 
constants by non-linear squares, the models were reparameterized as suggested in Draper 
& Smith by shifting the rate constant k 2 by the factor of 1000/T-3.333. The preexponential 
constant therefore, was centered at 300 K with the determined activation energy being 
1/1000 of the true activation energy. For K, 1/T was shifted to 1000/T-0.54, centering Ka 
at 1852 K. This resulted in faster converging least squares routines. A limit was also 
placed on the maximum value of D, the desorption activation energy, of 574 kJ/mole, 
preventing arithmetic overflows during the determination of the constants. 

The constants for Model 5 were determined using LOTUS 123 with a single parameter 
perturbation interactive approach. Each parameter was individually varied to determine 
the direction of the sum of squares of error (SSE) on log gamma. Then, depending upon 
which parameter reduced the SSE the most, that parameter was corrected and new 
parameters were found. The parameters for Model 6 were determined by choosing 
activation energy for the low temperature reaction coefficient, k 2 . The values used were 
640 cal/mole or 2.67 kJ/mole for nitrogen, and 1,000 cal/mole or 4.19 kJ/mole for oxygen 
atoms. These numbers are based on data for recombination of atoms on silica as 
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presented in Newman’s and Seward’s theses (1987, 1985). The remaining parameters in 
that model were determined by non-linear least squares. 

A pplication of Models 

Table 4 presents the overall summary. The models are compared to the recombination 
data for the recombination of nitrogen and oxygen atoms on RCG coated HRSI materials. 
The best fits were achieved by Model 1, as demonstrated by the low SSE in Table 4 for 
both cases. Willey’s Model 6 also closely resembles the data, giving a slightly lower SSE 
for nitrogen and slightly higher SSE for oxygen. Figures 6 through 12 and Figures 13 
through 19 (beginning on page 32) show a graphical representation for various models 
compared to the data used. Figures 6 through 12 illustrate recombination of nitrogen 
atoms on RCG coated HRSI surfaces and Figures 13 through 19 illustrate recombination 
of oxygen atoms on RCG coated HRSI surfaces. This arrangement portrays the relation 
of each model to the data for recombination of nitrogen and oxygen atoms. Figures 20 
through 23 are residual plots for the various models evaluated: Figure 20 for the 
recombination coefficient for nitrogen as a function of 1/T; Figure 21 for the nitrogen 
recombination coefficient versus partial pressure of nitrogen atoms; Figure 22 for the 
various recombination models for oxygen atoms as a function of 1/T; and Figure 23 for 
the various recombination models for oxygen atoms as a function of the partial pressure 
of oxygen atoms. These residual plots are helpful in determining the sufficiency and 
insufficiency of the various models used by noting trends and making corrections as 
needed. 

The results from Model 0 are shown in Figures 6 and 13 for nitrogen and oxygen 
recombination coefficients respectively. Model 0 predicts a straight line on an Arrhenius 
type plot with deviations occurring at higher temperatures. Figures 20-A for nitrogen and 
22-A for oxygen show the lack of fit at low ITT values where the actual data falls below 
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TABLE 4 


Residual Sum of Squares Summary for Various Recombination Coefficient Models 


Model # 

Name 

S.S.E. ; for 
Nitroeen 

S.S.E. for 
Oxveen 

0 

Arrhenius 

0.254 

1.357 

1 

Rideal-Eley 
without Concen- 
tration Terms 

0.138 

0.938 

2 

Rideal-Eley 

0.150 

1.044 

3 

Rideal-Eley 
with Temperature 
Correction Factor 

0.167 

1.044 

4 

Langmuir- 

Hinshelwood 

2.259 

1.019 

5 

Seward’s 

0.206 

1.134 

6 

Willey 

0.136 

1.021 


1 Sum of Squares of Errors on log gamma 
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that predicted. It is necessary to create models that would better predict recombination 
coefficients at higher temperatures. 

The results of Model 1 are shown in Figures 7 and 14 for nitrogen and oxygen atom 
recombination coefficients. These figures show a maximum in recombination coefficient 
as temperature increases; the maximum occurs at about 1600 K with a maximum predicted 
recombination coefficient of 0.013 for nitrogen and of 0.014 for oxygen. The residual plots 
confirm the suitability of this model. These plots are shown in Figures 20-B, 21-B, 22- 
B and 23-B. Figure 23-B shows a slight negative trend, indicating that the model is 
deficient in a concentration term. 

The results for Model 2 are shown in Figures 8 and 15 for nitrogen and oxygen atom 
recombination coefficients, demonstrating the influence of atom concentration (atom partial 
pressure) in terms of the Rideal-Eley mechanism. Figure 8 shows predicted recombination 
coefficients for partial pressure of nitrogen atoms at the surface of 100, 500 and 2,000 Pa. 
The recombination coefficient has a first-order functionality on partial pressure 
(concentration) at higher temperatures; at low temperatures, the three pressure predictions 
combine into a single curve on the right hand side of Figures 8 and 15 because the partial 
pressure functionality is zero-order. Thus the Rideal-Eley model predicts a zero-order 
concentration functionality at lower temperatures. 

Model 2 shows wider residual plots (Figures 20-C, 21-C, 22-C, 23-C) compared to Model 
1, and a negative trend in pressure appears for oxygen atoms (Figure 23-C). Nitrogen 
residuals also show a slight negative trend in terms of pressure (Figure 21-C). These 
residual plots, as a function of pressure, indicate that the model should be at least one 
order less in pressure. 

In Model 3, the addition of a temperature concentration term (2Mif/k 5 T) ;/2 does not 
enhance the fit of data. Figures 9 and 16 show very similar results to Figures 8 and 15, 
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however a slight concave downward appearance occurs in the curves for the low 
temperature region between 300 and 1000 K. This appearance occurs because of the T ;/2 
in the numerator of the Model 3. Again, the recombination coefficient is independent of 
concentration up to about 1500 K, after which it becomes first-order in concentration. The 
first-order influence can be seen in Figure 9 where the curves separate at higher 
temperatures for the partial pressures of nitrogen of 100, 500, and 2000 Pa. A further 
complication occurs in the determination of the parameters for this model; the upper 
boundary of the adsorption activation energy of 574 kJ/mole was reached for both oxygen 
and nitrogen. As further supporting data becomes available, the maximum can be 
increased. 

Model 4, the Langmuir-Hinshelwood Model, shows a poor fit for the nitrogen data (Figure 
10) and a good fit for the oxygen data (Figure 17). In determining the parameters for 
nitrogen atom recombination in Model 4, boundaries were reached in which an arithmetic 
overflow occurred during the computational steps. Therefore, the parameters chosen are 
at the boundary for Ka and for D. The model essentially looks like a negative first-order 
for nitrogen atoms throughout the whole region. In Figure 17, the recombination rate is 
higher for lower partial pressures for oxygen atoms with the curves at 100, 400 and 800 
Pa falling beneath each other sequentially. The data, which had different pressures 
depending upon the researcher, were well described by Model 4 considering the wide 
range of variables studied. These results suggest that measurements be taken in the 500 
to 1000 K range for oxygen recombination, varying the oxygen concentration at two or 
three levels. For instance, experiments could be run in the current arc jet facilities 
studying oxygen concentrations of 5%, 10%, and 21%, which would vary the concentration 
of atoms at the surface. The Langmuir-Hinshelwood model explains the observation of 
higher recombination rates at lower pressures; a maximum in the recombination coefficient 
is reached in this model as a function of concentration because the reaction requires two 
neighboring sites. This maximum occurs when the surface is about half covered with 
adsorbed atoms, when half of the sites are empty and half of the sites are covered. When 
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the surface is bare no reaction occurs, because nothing is on the surface. On the other 
extreme, at higher pressures, where the surface is covered, atoms which are near the 
surface cannot adsorb but must wait until two atoms combine on the surface and desorb; 
therefore, the maximum occurs when the surface is half covered. Although good data is 
unavailable to confirm this model, experiments need to be conducted at temperature 
regions between 500 and 1000 K with the partial pressure of oxygen purposely varied in 
order to determine an influence of recombination coefficient as a function of pressure. 
These experiments would provide the evidence required to support a Langmuir- 
Hinshelwood recombination model. 

The results for Model 5 are shown in Figures 11 and 18 for nitrogen and oxygen atom 
recombination coefficients respectively. In Figure 11, as temperature increases, the model 
is concave upward. Seward’s model achieves this by changing the parameters P and S, and 
a maximum gamma is reached at about 1600 K. At temperatures above 1,600 K, gamma 
decreases and the functionality is first-order in the partial pressure of atom concentration. 

At high temperatures, Model 5 is very similar to the Rideal-Eley mechanism. The residual 
plots for Model 5, Figures 20-F, 21-F, 22-F, and 23-F, show good scatter of the data points 
around zero. 

Model 6, the empirical model, is shown in Figures 12 and 19 for nitrogen and oxygen atom 
recombination coefficients respectively. This particular model has no concentration term, 
but does predict that the recombination coefficient will increase to a maximum as 
temperature increases and then decrease as temperature increases further. Model 6 also 
demonstrates concave upward behavior in the low temperature region. This is based on 
data from atom recombination on silica type surfaces which are very similar to RCG 
coated surfaces. These results suggest that the activation energy at low temperature is 
lower compared to activation energy at higher temperatures. Thus, the activation energy 
changes with temperature. The observed activation energy is the combination of two _ 
parallel reaction’s activation energies; one reaction going faster at low temperature and 


RJW/cv/014 


28 


the other speeding up at high temperature. The activation energy chosen for the low 
temperature reaction was based on values reported by Newman (1987) and Seward (1985) 
for nitrogen and oxygen atoms respectively. They reported the activation energy to be 640 
cal/mole for nitrogen atom recombination, and 1000 cal/mole for oxygen atom 
recombination on silica surfaces at low temperatures. The residual plots for Model 6, 
shown in Figures 20-G, 21-G, 22-G, and 23-G, demonstrate a good fit in all cases, except 
for a slight negative trend in terms of pressure for oxygen. 

Table 5 presents the parameter estimates for the various recombination models evaluated. 
Column 1 represents the model number which corresponds to the model number in Table 
3, columns 2 and 3 present the parameter name and its units. Rate constants and 
equilibrium constants are composed of two terms: a preexponential constant and an 
activation energy. Columns 4 and 5 are the parameter values which give the minimum 
sum of squares of errors for the comparison of each model to the data. The confidence 
limits for the parameters are not presented because of the poorly conditioned non-linear 
surfaces encountered with these models. 

In summary, a model independent of concentration describes the available data. The 
recombination model can be one of two versions, depending upon future results of low 
temperature recombination experiments: if data falls in the lower regions at 500 K 
(Figures 12 and 19) Model 6 should be used; if data falls upon a straight line between 300 
and 1000 K (Figures 7 and 14) Model 1 is more appropriate. If, in a further work, a 
concentration functionality is found important, then data should be compared to Models 
2 to 5 (Figures 8 to 11, Figures 15 to 18). 
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TABLE 5. 


Parameter Estimates for Various Recombination Models Evaluated 


Model # 

Parameter 

Name 

Parameter 

Units 

Value for 
Nitrogen 

Value for 
Oxvgen 

0 

^20 

Unitless 

0.019879 

0.017057 


Ea 

kJ/mole 

13.247 

12.585 

1 

^20 

Unitless 

0.021953 

0.020155 


Ea 

kJ/mole 

13.580 

13.235 


Ka 

Unitless 

1.0647E-07 

5.7376E-16 


D 

kJ/mole 

249.318 

520.880 

2 

^20 

unitless 

0.021660 

0.019763 


Ea 

kJ/mole 

13.545 

13.186 


Ka 

m J /atom 

1.2309E-33 

8.6919E-40 


D 

kJ/mole 

359.890 

574.340 

3 

k-20 

m/sec 

9.64640 

8.54460 


Ea 

kJ/mole 

15.971 

15.705 


Ka 

m J /atom 

9.1630E-35 

9.7166E-40 


D 

kJ/mole 

404.250 

574.340 
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TABLE 5. Continued 


Model # 

Parameter 

Name 

Parameter 

Units 

Value for 
Nitrogen 

Value fpr 
Oxygen 

4 

^20 

atom/m 2 /sec 

1.5283E+23 

8.1484E+22 


Ea 

kJ/mole 

14.209 

12.676 


Ka 

m J /atom 

1.9649E-19 

1.8112E-38 


D 

kJ/mole 

574.340 

574.340 

5 

D 

kJ/mole 

353.500 

339.000 


Ea 

kJ/mole 

4.593 

4.185 


Q 

atoms/m 2 

1.10E+20 

7.94E+19 



unitless 

2.25E-04 

1.79E-04 


P; 

1/K 

3.12E-03 

2.92E-03 


Max P 

unitless 

0.04 

0.10 



unitless 

3.16E-01 

5.01E-01 


Sy 

1/K 

-2.00E-03 

-2.00E-03 

6 

kio 

unitless 

3.9885E-02 

7.4805E-01 


Ea 2 

kJ/mole 

20.409 

58.861 


kj<j 

unitless 

2.4275E-04 

4.2845E-04 


Ea 2 

kJ/mole 

2.677 

4.189 


Ka 

unitless 

2.8098E-08 

1.4158E-14 


D 

kJ/mole 

262.921 

458.036 
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FIGURE 7 

Comparison of Model 1 with the Recombination Coefficients 
for Nitrogen Atoms on RCG Coated HRSI Materials 


Model 1 - Gamma N 
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FIGURE 8 




Comparison of Model 2 with the Recombination Coefficients 
for Nitrogen Atoms on RCG coated HRSI Materials 


Model 2 - Gamma N 
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FIGURE 10 

Comparison of Mode] 4 with the Recombination Coefficients 
for Nitrogen Atoms on RCG Coated HRSI Materials 


Model 4 - Gamma N 
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FIGURE 11 


Comparison of Model 5 with the Recombination Coefficients 
for Nitrogen Atoms on RCG Coated HRSI Materials 


Model 5 - Gamma N 



1000/T, i/K 
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FIGURE 12 


Comparison of Model 6 with the Recombination Coefficients 
for Nitrogen Atoms on RCG Coated HRSI Materials 


Model 6 - Gamma N 



FIGURE 13 

Comparison of Model 0 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 
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FIGURE 14 


Comparison of Model 1 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 
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FIGURE 15 


Comparison of Model 2 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 


Model 2 - Gaia 0 



1000/T, i/K 
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FIGURE 16 

Comparison of Model 3 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 


Model 3 - Gama 0 
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Comparison of Model 4 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 


Model 4 - Gamma 0 







FIGURE 18 

Comparison of Model 5 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 


Model 5 - Gamma 0 




FIGURE 19 

Comparison of Model 6 with the Recombination Coefficient 
for Oxygen Atoms on RCG Coated HRSI Materials 


Model 6 - Gamma 0 
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Residual Plots for Various Recombination Models 
for Nitrogen Atoms as a Function of 1/T 
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FIGURE 20 G 


Residual Plots for Various Recombination Models 
for Nitrogen Atoms as a Function of 1/T 
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FIGURE 21 A-C 

Residual Plots for Various Recombination Models 
for Nitrogen Atoms as a Function of Pressure 
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Residual Plots for Various Recombination Models 
for Nitrogen Atoms as a Function of Pressure 
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FIGURE 21 G 
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FIGURE 22 A-C 


Residual Plots for Various Recombination Models 
for Oxygen Atoms as a Function of 1/T 
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FIGURE 22 D-F 


Residual Plots for Various Recombination Models 
for Oxygen Atoms as a Function of 1/T 
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FIGURE 23 A-C 

Residual Plots for Various Recombination Models 
for Oxygen Atoms as a Function of Pressure 
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Applications of Models 1 and 6 to Heating Rate Determination for the Aeroassist Flight 
Experiment 

The application of recombination Models 1 and 6, for the aeroassist flight experiment, to 
heating rate design as a function of time for reentry trajectory was computed by Mr. Stan 
Bouslog of Lockheed Engineering and Sciences Company. Comparing the predicted 
heating rates of Kolodziej and Stewart’s (1987) catalytic models, which are discontinuous 
as a function of temperature, Models 1 and 6, which are continuous models, present 
slightly higher heating rates earily in the reentry pattern. However, Models 1 and 6 
predict much lower heating rates when compared to the fully catalytic wall (shown by the 
solid curve in both figures). In time, a brief transition occurs where Models 1 and 6 
predict a lower heating rate than those predicted by Kolodziej and Stewart due to their 
lower recombination coefficient. The peak heating rate in both cases occurs about 110 
seconds into the reentry trajectory. The maximum heating rate for Models 1 and 6 is 
about 35 BTU/ft 2 /sec, compared to the maximum heating rate of about 33 BTU/ft 2 /sec 
predicted by Kolodziej and Stewart. Thus, the goal of providing design engineers with a 
continuous recombination model for atom recombination on surfaces has been achieved. 
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Comparison of Stagnation Point Convective Heating Rates for 
Catalysis Model 1 Compared to Fully Catalytic and Kolodziej and Stewart’s Reaction 

Models 
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CONCLUSIONS 


Recombination models for oxygen and nitrogen atoms on surfaces can be derived from the 
Rideal-Eley mechanism which describes observed recombination data on RCG coated 
HRSI surfaces. The most appropriate model was dependent in temperature only and 
included an adsorption equilibrium term in the numerator and denominator. Normally, 
Rideal-Eley rate models for recombination coefficients are zero-order in concentration at 
low temperatures, and first-order at high temperatures; however, available data does not 
support inclusion of concentration terms into recombination coefficient models. Two 
recombination models were applied to stagnation point heating rate simulations and 
illustrate similar heating rates to those achieved by catalytic relations discontinuous in 
temperature. Parameter estimation in all of the models is difficult, requiring 
reparameterization and possible pre-determination of some constants. More experiments 
are recommended to verify these conclusions. 

RECOMMENDATIONS 

More experimental work is recommended in two temperature ranges, 300 to 1000 K and 
1500 to 1800 K, regarding atomic recombination rates on RCG coated HRSI materials. 
Additional experiments which vary the partial pressure of atoms at the surface are also 
necessary. In summary, depending upon the results, one of the four models presented 
above can be used. If additional experiments support no concentration dependency, then 
Model 1 or Model 6 should be chosen (Model 6 if the apparent activation energy 
increases in the temperature range from 300 to 1400 K, Model 1 if the constant activation 
energy occurs from 300 to 1400 K). If concentration dependence is found, then Model 
2 or Model 4 should be chosen (Model 2 if the recombination coefficient is found to be 
zero-order at low temperatures and first-order at high temperatures, Model 4 if the 
recombination coefficient is found to be negative first-order and transitions to first-order 
as temperature increases). 
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T.TST OF NOMENCLATURE 

Area normal to flux of atoms, m 2 
Average velocity of a gas molecule or atom, m/sec 
Concentration of surface sites in Seward Model, sites/m 2 
Equilibrium activation energy for adsorption, J/atom or kJ/g mole 
Activation energy for Rxn.-A-2, J/atom or kJ/mole 
Planck’s constant, 6.6256 E-34 N-M sec 
Equilibrium constant for Reactions A1 & B-l, m 5 /atom 
Arrhenius preexponential constant for Equation la, unitless 
Arrhenius preexponential constant for k. 70 Rxn. A-l & B-l, 1/sec 
Forward reaction rate constant in Rxns. A-l & B-l, m 3 /sec site 
Forward reaction rate constant in Rxns. A-2 & B-2 m 5 /sec site 
Forward rate constant used in Model 6, unitless 

Arrhenius preexponential constant for k 2 in Rxns. A-2 & B-2 m /sec site 

Arrhenius preexponential constant for k 3 , unitless 

Arrhenius preexponential constant for K in Rxns. A1 & B1 m J /atom 

Boltzmann Constant, 1.38 E-23 N-m/K atom 

First-order rate constant, m/sec 

Mass of an atom, kg/atom 

Nitrogen atoms 

Concentration of nitrogen atoms, atoms/m 5 7 

Surface impingement rate of atoms as given by kinetic theory, atoms/nrsec 
Oxygen atoms 

Concentration of oxygen atoms, atoms/m 5 
Concentration of adsorbed oxygen atoms, atoms/m 2 
Steric factor in Eqn. A- 19 = P 0 exp (P 2 T) 

Preexponential factor for P 
Exponential factor for P 
Partial pressure of gas species i, Pa 
Total pressure, Pa 

Universal Gas Constant 8,314 kJ/kg mole/K or 0.0008314 kJ/g mole K 
Concentration of reduced sites on the surface, atoms/m" 

Concentration of total sites on the surface, atoms/nr 
Preexponential constant in So 
Exponential factor in So 

Sticking coefficient in Eqn. A-19 = S 01 exp (S 7 T) 
time, s 

Temperature at surface, K 
Temperature at the tile surface, K 
delta, thermal desorption coefficient in Eqn. A-19 
gamma, recombination coefficient 
Stefan-Boltzmann constant 
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APPENDIX A 


Derivation nf Rideal-Elev Mechanistic Models 
The basic rate equations are: 

Adsorption of oxygen onto the surface: 

kj 

O + [S] ^ [OS] A-l 

k; 

Reaction of a surface oxygen with a gas phase oxygen: 

^2 

Og + [OS]— > O, + [S] A-2 

Where [S] is the concentration of "reduced" surface sites, sites/m 2 

and where [OS] is the concentration of "oxidized" surface sites, atoms/m 2 . 

Based on the law of mass action, the rate of 0 2 production per unit area is proportional 
to the concentration of oxygen in the gas phase times the concentration of oxide sites: 

d 0 2 /A„dt = k 2 [O] [OS] A- 3 

where [O] is in atoms/m 3 and k 2 has the units m 3 /sec atom. 

At steady state the rate of surface oxide, [OS], site formation has to equal zero: 

d[OS]/A*dt = 0 = k; [0][S]-k i [0S]-k 2 [0][0S] A-4 

If the number of total surface sites are constant then 

[S„]=[S] + [OS] a-5 

or 

[S] = [S„]-[OS] A-6 

Substitution of Equation A-6 into Equation A-4 leads to the following equation: 

d[OS]/A„ dt=0=k ; [O][S 0 ]-k / [O][OS]-k / [OS]-k 2 [O][OS] A-7 

[OS] can be factored and a solution found for [OS] in terms of rate constants, gas phase 
concentration, and catalyst site concentration results. 
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[OS] = kj[0][S (J ]/(k J [0]+k, ; +k 2 [0]) A-8 

Substitution of Equation A-8 into Equation A-3 leads to: 

d0 2 /A N dt=k 2 k 2 [0] 2 [S,]/(k ; [0]+k. ; +k 2 [0]) A- 9 

Gamma 0,)f 0 , is defined as the number of atoms reacting on the surface divided by the 
number of atoms which strike the surface: 


5f„=-dO/A„ dt/N c 

A-10 

By reaction stoichiometry: 


O+O ^ 0 2 

A-ll 

j 

or in terms or rates, 


-1/2 dO/dt = dO,/dt 

A-12 

or 


dO/dt = -2 dO,/dt 

A-13 

Substitution of Equation A-13 into Equation A-10 results in: 


to = 2 d0 2 /A„ dt/N 0 

A-14 

or 


to = 2 kjfclOl^SoKk.IOJ+kj+kjlOD/No 

A-15 

From the kinetic theory of gases 


N 0 =[0]c/4=[0][k/T/2mtf] ;/2 

A-16 

where 


c is the average velocity per particle = [8k s T/m’)f # ] ;/2 

A- 17 


and 


[O] is the number density of atoms, atoms/m 5 
Substitution of Equation A-16 into Equation A- 15 results in the following: 
^=2k i k 2 [0][S,][2m'7T / /k fl T] ;/2 /(k J [0]+k 2 +k 2 [0]) 


A-18 
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Discussion: Equation A-18 is analogous to recombination models derived by Seward (1985) 
for oxygen recombination and Newman (1987) for nitrogen recombination. 

Seward’s model is as follows: 

y 0 =2PN 0 So exp(-Ea/k 5 T)/(N 0 So +£+PN 0 exp(-Ea/k B T) A-19 

By factoring [O] in Equation A-18 and N 0 in Equation A-19 the equations can be 
compared. 


Y 0 =2k ; k,[So](2m /VTy /2 /(k ; +k,/[0]+k 2 ) A-20 


and 


'# 0 =2 P So exp (-Ea/k s T)/(So +®/N 0 + exp (-Ea/k fl T)) A-21 
Thus the adsorption term in Equation A-l is equivalent to So 

k, = So A-22 

likewise 

MO] = S/N 0 A-23 

and 

k, = P exp (-E/k/T) A-24 

In conventional kinetics, each rate constant, k„ can be described by an Arrhenius 
expression: 


k,.=k,. 0 exp (-E/k/T) A-25 

A comparison of Equation A-20 and A-21 can be made by placing the Arrhenius 
definitions on the left hand side of equations A-22, A-23 and A-24 and the definitions 
given in Seward (1985) and Newman (1987) for the terms on the right hand side. Thus, 

k i(3 exp (-E ; /k B T) = S 0] exp (-0.002 T) A-26 

The left hand side predicts that will increase with temperature while the right hand side 
predicts that will decrease slightly with temperature. This is t he major difference 
between the two models. 
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k 10 exp (-D/k 5 T)/[0]=C fl (k s T/h)exp(-D/k B T)/N 0 


A-27 


or 

^^“C^kjT/hyiOK^TCm®" 2 A-28 


or 

k 10 =C a (k B T2rrfJ7) ;/2 /h A-29 

Equal within a constant with a small influence of temperature 
and 


k 20 exp (-Ea/k fi T)= P exp (-Ea/k B T) 


A-30 


or 


k 2 o — P 

In Seward’s work P is a function of temperature P 0 exp (P 2 T) until P exceeds 0.1 where 
it is set equal to 0.1. Therefore, k 20 is within a constant factor of P at higher 
temperatures. 


Simplification of Equation A-18 

Often the temperature term [2m‘^7k fl T] ;/2 is assumed approximately constant in the 
temperature range of interest as compared to the exponential terms, and is combined with 
a preexponential constant in the Arrhenius equation. If the rate of reaction A-2 is rate 
limiting, k 2 [Oj is much lower compared to k ; and k 2 [O] in the denominator and the term 
is dropped. Dividing top and bottom by k. ; , one gets the classical Rideal-Eley mechanism 
with Equation A-2 as the rate limiting step. 

8 0 =2k 2 K[0]/(l+K[0]) A-31 

where K is the equilibrium constant for Reaction A-l. 
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Another assumption assumes that oxygen atom density is approximately constant near the 
surface and oxygen atom density can, therefore, be combined with the preexponential 
constants of the Arrhenius equation: 


'S a =2k 2 K/(l+K) 

A-32 

This is the simplest 4 parameter model possible. 


where k 2 = k 20 exp (-Ea/k B T) 

A-33 

and K = Ka exp (D/K fl T). 

A-34 


Note: In the report, R, the universal gas constant, is used in place of k B . 
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APPENDIX B 


Derivation of a Langmuir Hinshelwood Mechanistic Model 
The basic rate equations are: 

Adsorption of oxygen onto the surface 

k; 

o + [S]^=* [OS] B-l 

k.j 

Reaction of 2 surface oxygens to form molecular oxygen 

[OS] + [OS] — *> 0 2 + 2[S] B-2 

In this derivation Reaction, B-l is assumed in pseudo-equilibrium and Reaction B-2 is rate 
limiting, thus: 

d 0,/A„dt = k,[OS]- B-3 

and 

k 2 [0][S]=k.,[0S] B-4 

If the number of total surface sites are constant, then: 

[S 0 ]=[S] + [OS] B-5 

or 

[S]=[S„]-[OS] b-6 

Substitution of Equation B-6 into Equation B-4 leads to the following equation: 
k 2 [0][S„] = k_, [OS] + k, [0][0S] B-7 

or solving for [OS], 

[OS]= k^OHSJ/Ckj+k^Ol) b-8 
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Multiply top and bottom by k 1 and define K to be the equilibrium constant, Equation B- 
8 becomes: 


[OS] =K[O][S 0 ]/(l + K[OJ) B-9 

Substitution of Equation B-9 into Equation B-3 leads to: 

d0 2 /A^t=k 2 K 2 [0] 2 /(l+K[0]) 2 B-10 

From Appendix A above, 

2 d0 2 /dt=-dO/dt B-ll 

and 

N 0 =[0][k fi T/2m]j] 7/2 B-12 

thus, 

y o = 2 k,[2m^/k B T] //2 K 2 [0]/(l + K[0]) 2 B-13 

The rate constant, k 2 and equilibrium constant, K, follow Arrhenius equations. 

k 2 =k 20 exp (-Ea/k^T) B-14 

K=Ka exp (D/k 5 T) B-15 
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APPENDIX C 


Data Base Used in the Study 
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